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ABSTRACT 

I propose a scenario for the formation of double-ring systems, as observed in some planetary 
nebulae and in the outer rings of SN 1987A. In this scenario, two jets, one jet on each side of 
the equatorial plane, expand into a thin, dense shell. Such a flow is expected in binary systems 
where the mass-losing primary undergoes an impulsive mass loss episode which forms a thin, dense 
shell. I assume that a small fraction of that mass is accreted onto a companion which blows the 
jets. Each jet accelerates shell's material it hits sideways, forming a higher density ring. Using 
several simplifying assumptions, I derive an expression for the radius of the ring, which depends 
relatively weakly on the jets and impulsive mass loss episode properties. This shows that such a 
scenario is feasible. If there are several such impulsive mass loss episodes, more double rings will 
be formed, as observed in the planetary nebulae MyCn 18 (the Hourglass nebula). Because of the 
binary interaction and orbital motion, the double-ring system is displaced from the symmetry axis 
of the main nebula, as observed. 

Subject headings: stars: mass-loss — binaries: general — circumstellar mater — planetary 
nebulae: general — supernovae: individual (SN 1987A) 
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1. INTRODUCTION 

Two rings, one on each side of the central star (hereafter termed "double rings"), have been 
reported for several evolved systems. The most famous is the double-ring system (the outer rings) 
of SN 1987A (Panagia et al. 1996; Crotts, & Heathcote 2000; Maran et al. 2000, and references 
in these papers). Double-ring systems have been found in the planetary nebulae (PNe) Hubble 
12 (Welch et al. 1999) and He 2-113 (Sahai, Nyman, & Wootten 2000). The PN MyCn 18 (the 
Hourglass nebula; Sahai ct al. 1999) contains several pairs of rings, and it is not clear yet if their 
formation mechanism is the same as that of double-ring systems. 

Recently, Hony et al. (2001) reported the discovery of a double-ring structure around the 
central star of the r] Car nebula. However, higher resolution observations show that these rings are 
not real (Smith et al. 2001b). There are other arguments contradicting the claim by Hony et al. 
(2001). (i) The more prominent ring of rj Car (to the southwest) has a circular appearance, or it 
is even slightly elongated in the radial direction from the central star (20/xm map of Hony et al. 
2001). However, it is expected that circular rings should look like ellipses elongated in the azimuthal 
direction, as in SN 1987A (Crotts & Heathcote 2000) and the several rings in the Hourglass nebula 
(MyCn 18; Sahai ct al. 1999). (ii) The two rings of r] Car are projected on the southeast lobe, 
which is tilted toward us (Davidson et al. 2001). There are several fainter 'blobs' and 'fingers' 
seen on the other lobe, which are generally symmetric with the two rings on the southeast lobe. 
It is quite possible that if the northwest lobe had the same inclination as the northeast lobe, the 
blobs and filaments would be brighter. It is possible, therefore, that the 'rings' are a result of some 
instabilities, or stochastic mass ejection on the northwest lobe, (iii) The long-slit spectroscopy used 
by Davidson et al. (2001) to infer velocities on the Homunculus covers a portion of the southwest 
ring. No peculiar velocities are observed there. This adds support to the idea that the 'rings' are 
located on the lobe, (iv) Another possibility is that the rings are in the orbital plane; the bipolar 
nebula around the luminous blue variable G25.5-I-0.2 is similar to the Homunculus of rj Car, and has 
a double-nested equatorial ring system (Clark, Steele, & Langer 2000). For these reasons and the 
evidence from the new observations of Smith et al. (2001), I will no longer refer to the previously 
claimed double-ring system in r] Car. 

Some interesting properties of the double-ring systems are; (1) The outer rings of SN 1987 are 
much more prominent than any of the other double-ring systems, with their density being much 
higher than the density of their surroundings (Burrows et al. 1995). In the other systems, the 
medium around the rings can be seen, with a much lower density contrast than that in SN 1987A. 
For that reason, and because of the presence of the inner ring in the equatorial plane of SN 1987A, 
it is not clear whether the formation mechanism of the outer rings of SN 1987A is the same as that 
of double rings in PNe. (2) While in some systems (SN 1987; Hb 12) the line connecting the two 
rings' centers (hereafter double-ring axis) is more or less along the major (symmetry) axis of the 
nebula, in He 2-113 this line is highly inclined to the longer, presumably the symmetry, axis. (3) 
In all systems the double-ring system is displaced from the central star, i.e., the double-ring axis 
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misses the central star (SN 1987A; He 2-113) and/or the intensity around the rings is not constant 
(e.g., He 2-113; Hb 12). (4) The rings are detected via different emissions bands, e.g., the double 
rings of HB 12 are seen in [Fe II] (Welch et al. 1999), and those of He 2-113 in Ha. 

Many models were proposed for the formation of the oTitcr rings of SN 1987A, but most of 
them were found to have some problems. Models based on interacting winds (e.g., Blondin, & 
Lundqvist 1993) were criticized by Burrows ct al. (1995) and Meyer (1997), while the ionization 
model of Meyer was criticized by Soker (1999). Crotts & Heathcote (2000), discuss these models, as 
well as the one based on a combination of interacting winds and ionization (Chevalier &; Dwarkadas 
1995), and the model based on an equatorial pulsed mass loss (Soker 1999), and find that no model 
is capable of explaining all properties of the outer rings of SN 1987A. A recent model proposed 
by Tanaka & Washimi (2002), which is based on magnetic tension, has too many assumptions 
concerning the mass loss geometry and magnetic field, and it is not clear whether it reproduces the 
properties of the rings. 

In the present paper, I explore a different model for the formation of a double-ring system, 
which can also explain several pairs of rings. I examine a short pulse of mass loss from the primary 
mass-losing star, which forms a spherical, dense, thin shell. I assume, then, that a jet hits and 
interacts with the shell (one jet on each side of the equatorial plane). This type of flow is expected 
to occur if part of the mass lost to the shell is accreted by a companion, which then blows the two 
jets. The central star of He 2-113 is a Wolf-Rayet type star (a [WC] star); it is possible that an 
impulsive mass loss episode was connected with the transition to a [WC] star. Jets, one on each 
side of the equatorial plane (if they are not well-collimated, they are called CFW for "collimated 
fast wind"), which last for a long time are thought to lead to the formation of bipolar lobes (Morris 
1987; Soker & Rappaport 2000; Soker 2002). Accreting close binary companion can form bipolar 
PNe (Soker &; Rappaport 2000), i.e., those with two polar lobes and an equatorial waist between 
them, while wider companions may form elliptical PNe (Soker 2001), i.e., PNe with a large-scale 
elliptical shape, but which can possess small lobes or jet-like features. There is no sharp boundary 
between the two groups; the approximate orbital separation between the binary systems that form 
bipolar, and those which form elliptical, PNe depends on the mass and nature of the accretor, and 
on the properties of the wind (particularly its velocity), blown by the mass losing star (Soker 2001). 

The model is discussed and shown to be feasible in §2. Disks can precess, explaining the 
inclination of the double-ring axis with respect to other structures in the nebula. This, and other 
possible explanations are discussed in §3. In §3, I also discuss the displacement of the double-ring 
axis from the central star. I summarize the main results in §4. 

2. THE PROPOSED MODEL 

I treat a jet expanding into a thin shell, and accelerating the material in the shell, both 
forward and to the sides. I assume that the jet is formed by an accretion disk (Livio 2000), which is 
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formed via mass transfer from an evolved star into a compact companion. The scenario can work in 
principle in cases where the jet is formed via a single star mechanism, e.g., Garcia-Segura, & Lopez 
(2000). The binary system destroys the axisymmetrical nature of the flow. However, to facilitate 
an analytical treatment, I assume axisymmetry. I return to this point in §3.2. It is assumed that 
prior to and after the impulsive mass loss episode, the mass loss rate from the primary was too low 
for the formation of jets (Soker k. Rappaport 2000). The impulsive mass loss rate into the shell is 
Mp, and that into one jet is Mj\ their expansion velocities are Vp and Vj, respectively. The opening 
angle of the jet, from its symmetry axis to its edge, is a, i.e., the jet's head covers a solid angle of 
27r(l — cos a) ~ 7ra^, where I assume a ^ 1. The densities in the shell and jet are 

Pp = 1 — 2~' Pj - — r^~' 

respectively, and r is the distance from the center. The velocity of the jet head, Vh, is determined 
by equating the slow wind pressure on the jet head with that of the jet material. Since thermal 
pressure can be neglected in the pre-shockcd media, the expression is Pj{vj — v^)'^ = Pp{vh — Vp)^. 
For the mechanism proposed here to work, the jet needs to deposit its energy in the shell. I 
therefore assume that the jet substantially slows down in the shell, such that Vh <C Vj\ hence Vh is 
approximately given by 



Vh ^ Vp 



X 1/2' 
Vj Pj' 



1 + 

Vp \Pp 




The second equality comes from the requirement that the jet be energetic enough to form the rings. 
This can be seen by substituting typical values, 

^"^■^lyio^J [omMp) ^io°J ■ 

This scaling may result, for example, from a companion that accretes ~ 10% of the mass lost in the 
pulse and blows 10% of that into each jet, which has an opening angle of ~ 10°. In progenitors of 
PNe the shell velocity can be ~ 10 km s~^, while the jet velocity from a main sequence companion 
can be ~ 200 — 500 km s~^, and an order of magnitude faster from a white dwarf companion. 

Because the interaction takes place close to the central star where densities are high so cooling 
is fast (Soker 2002), and the shell is thin, no long lived hot bubble is formed, and the acceleration 
of the shell material by the pressure formed by the interaction lasts for a short time, tacc- I assume 
that as the jet interacts with the shell, a region with pressure equal to the ram pressure of the 
jet is formed. With the assumption <C Vj this ram pressure is P ~ PjVj. The thickness of the 
shell formed by the impulsive mass loss event is dp = Vptp, where tp is the duration of the pulse. I 
approximate the segment of the shell in the jet's vicinity as a plane, with the z coordinate along 
the symmetry axis, i.e., perpendicular to the plane of the shell, and with the y coordinate as the 
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distance from the symmetry axis in this plane. The mass of the shell per unit length along the z 
direction and in a circular cross section between the symmetry axis and a distance y from it, is 

m{y) = = Pp- (4) 

Let the shell material be accelerated to a transverse velocity vq as it leaves the interaction region 
at a distance yo = r sin a from the symmetry axis. I assume that the pressure accelerates the shell 
material from the symmetry axis {y = 0) to the jet head's radius yo, with a force per unit length 
of dF/dz = 2'iTyP. Conservation of momentum gives 



^^0^ Sc^(Pl] Vj5^vh5. (5) 



miyo) 



Pp) 



The transverse acceleration does not occur along the entire jet's cross section, hence it is less 
efficient than assumed. This is taken into account in the factor 5 < 1. 

The formation of a dense, thin envelope around jets propagating in a continuous medium 
is well established both numerically and analytically (e.g., Blondin, Fryxell, & Konigl 1990; de 
Gouveia dal Pino & Benz 1993; Cerqueira, de Gouveia dal Pino, & Herant 1997; Lee et al. 2001; 
Ostriker et al. 2001). In the present case, the medium is not continuous, hence a ring, instead of a 
thin envelope, is formed around the jet. Ostriker et al. (2001) derive an analytical expression for 
the shape of the thin envelope (which they call a "shell" ; not to be confused with the shell used 
here). Using their equation (22) close to the jet, for a short time after a parcel of gas leaves the 
interaction region, I derive the expression vq = {j3cs/vh)Vh, where vq and Vh are the notation used 
in the present paper, /? ~ 4 is a factor used by Ostriker et al. (2001) in their expression for the 
momentum flux, and Cg is their sound speed in the interaction region. Their ratio ^Cg/vh plays the 
role 5 plays here. They find the range of this ratio is 0.2 — 0.6, therefore I take 5 ~ 0.2 — 0.6, as 
well. Using equations (5) and (3), then, gives 



^0 



V400 kms-V VlOkms-V lo.OlMpi VlOV \OAJ 



The condition for equations (5) and (6) to hold is that the acceleration time assumed here, 
^acc ~ Vq/vq-, is shorter than the time it takes the jet to cross the shell, dp/vh^ which gives the 
condition yo ^ dp. In the other limit, yo ^ dp, the acceleration time is ^acc ^ dp I V}i , from which 
the velocity is derived 

~ Vh—^ for yo ^ dp. (7) 

However, in this case we also expect 5 to have a lower value. Hence, the case of a too thin shell and 
a too wide jet is not efficient for the proposed mechanism. The situation is more complicated. At 
large distances from the star (large r) the jet head cross section is larger, reducing the efficiency. 
On the other hand, at larger distances the radiative cooling is longer, making the acceleration more 
efficient (Soker 2002). 
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For the compression of rings to higher densities than their environment, the flow should start 
as supersonic vq > c, where c is the sound speed in the shell. Before ionization starts, the material 
blown by the primary evolved star is cool, T < 1000 K, and mainly atomic, or even molecular, 
with a sound speed of c < 4 km s~^. The interaction between the transversely expanding ring and 
the mass it sweeps from the cold shell is radiative, and the velocity is determined by momentum 
conservation, i.e., movo ~ m{y)v{y), where mo is the mass per unit length pushed by the jet to 
the sides. Not all of the shell mass that is hit by the jet head is pushed to the side; some is 
pushed forward, giving mo < m{yo). This is parametrized by mo = rym(yo)- For example, if the 
mass expelled to the sides comes from the region 0.8yo < y ^ Ho, then r] ~ 0.4. For y ^ yo the 
mass per unit length as given by equation (4) can be used. This gives v{y) = V'^^oiy /yo)~'^ ■ Using 
yo = rsina :^ ra gives the distance from the symmetry axis where compression stops as 

(8) 

r V c / 

With the aid of equation (6) this can be scaled as 

^ 3 ( \ ( ^ ( " ^ ( ] (9) 

r ■ UoOkms-V VlO km s"!/ U km s"!/ lo.OlMpi Vl0°/ \0A J \0AJ ' ^' 

The following should be noted in regard to the last equation. 

(1) This expression was derived under the assumptions that the shell segment from which the ring 
is formed is a plane, and the radial expansion (from the central star) of the material in the shell 
was neglected. Hence, this already approximate expression becomes very crude for large values of 

y/r. 

(2) The derived size of the ring for reasonable parameters in PNe is on the order of the observed 
size of the rings, with a relatively weak dependence on the different parameters. The same scaling 
may hold for the progenitor of SN 1987A during the formation of the outer rings, because the star 
was a red super giant, and the companion required by other arguments is a main sequence star in 
the mass range of 0.5 < M2 < 3Mq (Chevalier & Soker 1989; Podsiadlowski, Joss, & Rappaport 
1990; Collins et al. 1999). The outer rings were formed just before the companion entered the 
common envelope phase with the progenitor of SN 1987A. 

(3) The ring becomes larger for a larger opening angle a of the jet. However, a can't be too large 
because the conditions Vo > c and ra < dp {dp is the shell width) should be met. Therefore, the 
present scenario can't account for very wide rings close to the equatorial plane. Two such rings are 
observed in the massive eclipsing binary system RY Scuti (Smith, Gehrz, & Goss 2001a; Smith et 
al. 1999), and probably were formed by the equatorial flow (Smith et al. 1999). 

(4) The total interaction time is equal to the time it takes the jet to cross the shell dp/vh = 
tp{vp/vh) <^ tp {tp is the duration of the impulsive mass loss event that forms the shell). It turns 
out that the jet can be blown in a very short pulse. 

(5) As the jet breaks out of the outer side of the shell it pushes the dense shell's material to the 
sides and exterior to the shell (like a pencil punching a hole in paper) . This may be the explanation 
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for fainter gas extending radially from the rings in He 2-113 (see fig. lb of Sahai et al. 2000, 
particularly the eastern ring). We note that ionization fronts also tend to form denser tails behind 

dense material (see below). 

(6) Later evolution is not considered here. In particular, the central stars in all systems considered 
here become hotter, ionize the nebulae around them, and blow fast winds. Ionization proceeds faster 
in the low density surroundings of the ring, heating it to higher temperatures, hence increasing the 
surrounding pressure. The higher pressure may further compresses the ring and material in its 
shadow (e.g.. Canto et al. 1998; Pavlakis et al. 2001), increasing the density contrast. Later, the 
ionization heats the gas in the ring, which starts to expand, thus reducing the density contrast 
between the ring and its environment. The fast wind accelerates the lower density environment to 
high velocities, while the dense ring lags behind. This increases the density contrast. 

3. CHANGING THE SYMMETRY AXIS 

3.1. Changing the Direction of the Jets 

I consider two mechanisms for the inclination of the line joining the centers of the two rings to 
the symmetry axis of the main nebula, e.g., He 2-113 (Sahai et al. 2000). In the first, the accretion 
disk is precessing. The short life of the jets implies that they can be blown while the accretion disk 
has not get precessed very much. In the second mechanism, the impulsive mass loss of the shell 
was very clumpy. The rings in He 2-113 are also very clumpy (Sahai et al. 2000, fig lb). Because of 
the clumpy nature of the accreted slow wind, more angular momentum may be accreted from one 
side of the equatorial plane, and the accretion disk (if formed) will be inclined to the equatorial 
plane. As was found in the previous section, the mass loss rate into each of the two jets can be 
~ 0.01 of that in the impulsive mass loss event from the primary. Because a fraction of ~ 0.1 — 0.2 
from the accreted mass is blown into the two jets, only a fraction of < 0.1 from the mass lost by 
the primary needs to be accreted. This part can be clumpy. 

For accretion from a wind, the net specific angular momentum of the material entering the 
Bondi-Hoyle accretion cylinder with radius Ra = 2GM2/vr, i.e., the material having impact pa- 
rameter b < Ra, is Jbh = 0.5(27r/Po)-Ro (Wang 1981), where Pg is the orbital period, and Vr is 
the relative velocity between the wind and the companion of mass M2. Livio et al. (1986; see also 
Ruffert 1999) find that the actual accreted specific angular momentum for high Mach number flows 
is ja = VJEH, where r/ ~ 0.1 and r/ 0.3 for isothermal and adiabatic flows, respectively. Taking 
a circular orbit for simplicity, the relative orbital velocity of the two stars is Vorh = 27ra/Po) and 
taking 77 = 0.2 gives ja — O.lt'orb-^a/'^- '^^^ angular momentum is perpendicular to the equatorial 
plane, and for accretion from a homogeneous wind, the accretion disk around the companion will 
be in the equatorial plane. To significantly tilt the plane of the accretion disk, I consider a clumpy 
(inhomogeneous) wind, where an extra mass AM is accreted from one side of the plane. Let the 
wind speed be Vp as before, and assume that near the companion the extra mass has an impact 
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parameter of ~ 0.5Ra above the plane. The specific angular momentum added by this extra mass, 



where Mace is the total mass accreted to the disk. It is assumed that the companion accretes only 
0.1 of the mass lost by the primary. This implies that a > 0.5Ra- Also, for relevant PNe I find 
0.3 < Vp/vorh ^ 3. By using these values in equation (10) I conclude that the clumpy impulsive 
mass loss, with ^ 10% more mass blown on one side of the orbital plane in the companion direction, 
can significantly tilt the accretion disk and the jets blown by it. This requires that the impulsive 
mass loss episode, or at least the duration of blowing the jets, be much shorter than the orbital 
period. 



Several processes can lead to a departure from axisymmetry of nebulae formed by binary 
interaction (see §1 of Soker & Rappaport 2001). Here I briefly mention several effects, of which all, 
or a subset, of them can be significant in any particular case. The displacement of the rings can 
manifest itself in the rings' symmetry axis not crossing the central binary system, and/or one side 
of each ring (but the same side for both rings) being brighter or larger. 

Each of the two stars have an orbital velocity around the center of mass. If the impulsive 
mass loss episode is shorter, or not much longer, than the orbital period, then the mass leaves the 
system with a boost in the direction toward which the star moves at that period. This effect can 
influence the shell or the jet. The magnitude of the effect is approximately the ratio of the orbital 
to the expansion velocities. The shell is likely to acquire a larger departure from axisymmetry since 
its expansion velocity is much lower than that of the jets. 

Another effect is due to the companion being displaced from the center of the shell; the initial 
shell center is at the location of the primary mass losing star (at the moment it blows the shell), 
while the jets are blown somewhat later, and at the location of the binary companion. If the 

interaction between the jets and the shell occurs at a distance r, and the orbital separation is a, 
the magnitude of this effect is ~ a/r. Yet, in another process, the jets blown by the companion 
can be bent by the primary's wind (Soker & Rappaport 2000). 



Jclump) IS 




(10) 



3.2. Displacement from Axisymmetry 



I conclude that displacement from axisymmetry is a natural outcome of the scenario proposed 
in the present paper, and should occur in most cases. 
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4. SUMMARY 

A scenario for the formation of double-ring systems was analytically examined under several 
simplifying assumptions. In this scenario two jets, one jet on each side of the equatorial plane, 
expand into a thin, dense shell. Such a flow is expected in binary systems where the mass losing 
primary undergoes an impulsive mass loss episode which forms a thin, dense, and more or less 
spherical, shell, and a fraction ~ 0.01 — 0.1 of that mass is accreted onto a companion, such that an 
accretion disk is formed, and the companion blows the jets. To facilitate the derivation of a simple 
analytical relation that can be easily used, I approximated the segment of the spherical shell in and 
near the interaction region as a plane. I neglected the radial expansion of the shell, and assumed 
that a pressure equal to the ram pressure of the jet accelerates the shell material that is impacted 
by the jet. Material is pushed to the sides, emerging from the interaction region with transverse 
supersonic speed, vq > c. Some shell material is pushed forward (in the radial direction) by the jet. 
This reduces the efficiency of the acceleration, and lowers the mass of the shell that is hit by the 
jet and which flows to the sides. The reduced acceleration efficiency is represented by a parameter 
6 < 1, and the lower mass by a factor 77 < 1. 

The derived expression for the radius of the ring, y, is given in equation (9), which is scaled 
with typical winds' properties expected in progenitors of PNe. Equation (9) shows that the size of 
the ring depends relatively weakly on the winds' properties and the different physical variables and 
factors. Hence, I argue, the simple analytical treatment shows that such a scenario is feasible. The 
same scenario can account for the double-ring system in r/ Car (Hony et al. 2001), if real (however, 
sec §1). The double rings in SN 1987A (the outer rings) are very prominent. The process studied 
here by itself can't account for their formation. However, ionizing radiation emitted by, and a fast 
wind blown by, the progenitor of SN 1987A as it turned to a blue super giant, could account for the 
high density contrast between the rings and their surroundings. Three dimensional gas-dynamical 
numerical simulations are required to verify the proposed scenario. 

This scenario can account for several other properties of double-ring systems. If there are 
several such impulsive mass loss episodes, more double rings will be formed. Hence, this scenario 
may account for the rings in the PN MyCn 18 (the Hourglass; Sahai et al. 1999). The inclination of 
the symmetry axis of the double rings to the symmetry axis of the main nebula, which is observed in 
a few systems, can be explained by either a precessing accretion disk, or a clumpy (inhomogeneous) 
wind; accretion from a clumpy wind may form an accretion disk inclined to the equatorial plane 
(§3.1). The displacement of the double-ring systems from axisymmetry is observed in most cases. 
The displacement manifests itself in one side of the double-ring system being brighter or larger, or 
by their symmetry axis being displaced from the central star. This is expected in the presently 
proposed scenario. Because the impulsive mass loss episode is shorter, or not much longer, than 
the orbital period, the shell's center and jets' origin are in different locations. Another effect is the 
orbital motion of the two stars, which results in a net velocity of the shell and jets centers of mass 
relative to the center of mass of the binary system. The wind from the primary can bend the jets 
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blown by the secondary, in a process which also causes departure from axisymmetry. 
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